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Abstract —Radar imaging systems transmit modulated wide¬ 
band waveform to achieve high range resolution resulting in high 
sampling rates at the receiver proportional to the bandwidth of 
the transmit waveform. Analog processing techniques can be used 
on receive to reduce the number of measurements to N, the num¬ 
ber of potential delay bins. If the scene interrogated by the radar 
is assumed to be sparse consisting of K point targets, results from 
compressive sensing suggest that number of measurements can 
be further reduced to scale with KlogN for stable recovery of 
a sparse scene from measurements with additive noise. While 
unstructured random projectors guarantee successful recovery 
under sparsity constraints, they cannot be implemented in the 
radar hardware in practice. Recently, structured random Toeplitz 
and Circulant matrices that result from using stochastic wave¬ 
forms in time delay estimation setting have been shown to yield 
recovery guarantees similar to unstructured sensing matrices. 
However, the corresponding transmitter and receiver structures 
have high complexity and large storage requirements. In this 
paper, we propose an alternative low complexity compressive 
wideband radar sensor which combines multitone signal chirp 
waveform on transmit with a receiver that utilizes an analog 
mixer followed with a uniform sub-Nyquist sampling stage. 
We derive the recovery guarantees for the resulting structured 
measurement matrix and sufficient conditions for the number 
of tones. The only random component of our design is the 
sparse tone spectrum implementable efficiently in hardware. Our 
analytical and empirical results show that the performance of our 
scheme is in par with unstructured random sensing matrices and 
structured Toeplitz and Circulant matrices with random entries. 

Index Terms —Compressive sensing, mutual coherence. Struc¬ 
tured measurement matrix. Linear Frequency modulated wave¬ 
form, Radar. 


1. Introduction 

Radar imaging systems acquire information about the scene 
of interest by transmitting pulsed waveforms and analyze the 
received backscatter energy to form an estimate of the range 
and amplitude of the reflectors in the scene. These range 
profiles from multiple pulses and/or multiple antenna elements 
can be processed jointly to solve a multitude of inference 
tasks including detection, tracking and classification Q. In 
this paper, we focus on the problem of estimation of range and 
amplitude of refiectors in the scene using a single modulated 
wideband pulse (j)(t) of bandwidth B. The resolution of the 
echo imaging system is directly proportional to the bandwidth 
of the transmitted signal. Assuming the support of the observed 
delays are known to lie on an interval T (termed as range swath 
in radar literature), then the unknown range profile can be 
discretized into N = BT delay bins. The signal model at the 
receiver can be written as y(t) = + n{t), 


where A = l/B, n{t) is the receiver noise and Xn denotes 
the complex scattering coefficients associated with the n'th 
delay bin. Commonly the received signal is matched filtered 
with a copy of the transmitted pulse to detect refiectors in 
range and estimate their complex amplitude of backscattered 
energy. Direct digital implementation of the matched filtering 
step requires quadrature sampling of the received signal for 
the pulse duration with sampling rate matching the bandwidth 
of the transmit waveform. Alternatively, the matched filtering 
can be implemented in the analog domain where the number 
of samples are reduced to N to cover the delay support. 
However, the match filter output still requires to be sampled at 
the Nyquist rate corresponding to the system bandwidth, ren¬ 
dering digital and analog matched filter receivers for arbitrary 
waveforms impractical for bandwidths exceeding gigahertz for 
current analog to digital converter (ADC) technology. 

If a linear frequency modulated waveform (LFM) (j)(t) = 
is used on transmit, the matched filtering can be ap¬ 
proximately implemented through mixing the received signal 
with a reference LFM waveform and low pass filtering the 
mixer output. At the mixer output each copy of the waveform 
delayed by A appears as a tone whose frequency is given by 
This pre-processing step is termed stretch processing |[^, 
and can result in substantially reduced sampling rate for 
the ADC used in the mixer output if the delay support T is 
smaller than the pulse length. Specifically, received signal at 
stretch processor’s output can be written as: 

N 

= (1) 

n=l 

In essence, stretch processing converts range profile estimation 
problem into frequency spectrum estimation problem with 
Nyquist rate samples in time obtained after analog processing. 
If the scene can be assumed to have fewer targets K than 
the number of delay bins N, well-known results 0,0 from 
the Compressive sensing (CS) shows successful reconstruction 
with sub-Nyquist samples is possible with the number of 
measurements M scaling with KlogN, if appropriate mea¬ 
surement operators can be implemented. Furthermore, there 
are numerous tractable algorithms with provable performance 
that are either based on convex relaxation 0 0 or greedy 
methods 0 to solve the reconstruction problem. Moti¬ 
vated by these advances, compressed sensing techniques have 
been applied to a variety of problems in Radar | [T0| : range 
profile estimation m single pulse systems for range-doppler 
estimation in (13 single pulse multiple transmit and receive 
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system for range-doppler and azimuth estimation and target 
detection in (H), fill, r emote sensing in |T5| , direction of 
arrival estimation in 116|. CS based radar sensors based on 
pure random waveforms (II)> Xampling framework and 
Random Modulator Pre-Integrator (RMPI) framework |T9| 
using the receiver structure from p0| have also been im¬ 
plemented. A common theme in most of the CS literature 
has been randomization as it leads to measurement matrices 
that have provable recovery guarantees. Implementation of the 
randomness into compressive radar systems has proven to be 
a challenging task in practice; uncorrelated random signals 
with high peak-to-average power ratio is mismatched to the 
nonlinear power amplifiers used in radar systems and the 
system bandwidth (and as a result range resolution) is limited 
as digital to analog converters (DAC) have to be employed for 
generation of precise random signals for transmit and receive 
mixing. 

The LFM pulse model in Q provides an alternative strategy 
as it converts the range estimation problem into an equiva¬ 
lent sparse frequency spectrum estimation problem. Uniform 
subsampling in this setting has poor perforamnce 0- Non- 
uniform random sub-sampling can be used to obtain mea¬ 
surements with low mutual coherency 1^ , p2| . However, 
non-uniform sampling with commercially available ADCs still 
requires it to be rated at the Nyquist rate to accommodate 
closely spaced samples. We propose to use low speed uni¬ 
form sub-sampling using a high analog bandwidth ADC in 
the sparse frequency spectrum estimation setting and push 
randomness to transmit signal structure to obtain compressive 
measurements. ADCs whose analog bandwidth exceeds their 
maximum sampling rate by several factors is readily available 
and used routinely in pass-band sampling. This compressive 
radar structure proposed in uses linear combination of 
LFM waveform at the transmitter with randomly selected 
center frequencies, while maintaining the simple standard 
stretch processing receiver structure. The output of the stretch 
processor receiver is given by 


N Nc 

y{t) = X X 

n=l k=l 




where fn.k is a predetermined known complex phase, Nc is 
the number of tones modulating the LFM waveform. We ob¬ 
serve that under the proposed compressive sensor design each 
delayed copy of the transmitted waveform is mapped to multi- 
tone spectra with known structure. As shown in this paper this 
known multi-tone frequency structure enables recovery from 
aliased time samples with provable guarantees. These results 
complement previous work which has shown good empirical 
performance in simulation and measurements p4| . 

The contributions put forth in this sequel is that we present 
the theoretical analysis for this multi-frequency LFM system. 
We show that the system with a relatively small number 
of LFM waveform has performance guarantees similar to a 
matrix with independent random entries for a sufficiently large 
number of tones modulating the LFM waveform. We also 
present a numerical analysis comparing our system with other 
measurement schemes. 


Notation and Preliminaries 

We denote a vector in A/'-dimensional complex domain as 
X G c^. ||x||q is called as norm, which is given as the 
number of non-zero elements in a vector. Clearly, this is not a 
valid norm but is used in formulating the fundamental problem 
in compressed sensing. We denote ||x||^ = \xi\ as the 
norm. The Euclidean or ^2 norm is given as ||x ||2 = 

\j^i • We denote a matrix as A G A* G 

as conjugate transpose of A, and I as the identity matrix of 
dimensions dependent on the context. The spectral or operator 
norm of the matrix is given as ||A||^^ = crmax(A) the largest 
singular value of the matrix. The Frobenius norm of a matrix is 

given as ||A||^ = Another important quantity 

of interest is the mutual coherence /i (A), which is a measure 
of the correlation between the columns of matrix A. The 
mutual coherence is given as /i (A) = max^^j ||a^"||a^|| . 
where Ai is a column of matrix A. Another fundamental 
property for the measurement matrices is called as Restricted 
Isometry Property (RIP). A measurement matrix A G 
is said to satisfy RIP of order AT, if for any K-sparse vector 
xgC^ 

(l-,5)||x||2<||Aa:||2<(l + <5)||x||^ 
equivalently, Sk = max ||ApAr — I|| < (2) 

r 

card{V)<K 

where T is an index set that selects the columns of A, and 
cardiT) refers to the number of elements in the set, ^ < 1 
and Ar is the restriction of A having columns indexed by F. 
We denote the expectation operator as E(.). The circularly- 
symmetric complex Gaussian distribution with mean /i and 
variance is denoted as CM (/i, a^). 

A. Relation to other works 

Related results in literature on recovery guarantees for 
compressed radar sensing can be broadly categorized in two 
categories: Results relating to signal reconstruction establish 
uniform recovery guarantees for successful reconstruction of 
all K-sparse signals, whereas results on support recovery is 
concerned with the detection of non-zero locations of a K- 
sparse signal that is assumed to be sampled from a generic 
statistical model, such as uniformly sampling from all possible 
subsets of size K 

In this paper we show that randomly sampled K-sparse sig¬ 
nals can be recovered with high probability using LASSO for 
the structured measurement matrix of the proposed ompressive 
radar sensor sensing scheme. Next, we show that the estimates 
of mutual coherence and column norms we obtain can be used 
to provide uniform recovery guarantee following a standard 
argument. 

Table summarizes related results for support recovery 
for different measurement matrices. The upper bound on the 
sparsity level that guarantees successful support recovery for 
our scheme has an additional log (A + M) penalty compared 
to unstructured Gaussian matrix as shown by Candes and 
Plan in p6| and block Toeplitz matrices with entries sampled 
from Rademacher distribution as shown by Bajwa in We 
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Recovery Guarantees from noisy measurements with component-wise noise variance cr^ 

Matrix Type of size M x N 

Mutual Coherence 

Sparsity for successful recovery 

Spectral Norm 

Minimum signal 

strength 

Reference 

Random matrix with (NM) 
independent random entries 


°(n&) 

yi+i 

O (<tV2 log N) 

25 

- 

27 

Toeplitz block matrix with 
(N M) random entries 

o(yw) 

o(n*) 


0{aV2logN) 


28 


LFM waveform modulated 
with Nc ‘C N randomly 
selected tones 


of ^ \ 

^ \ \ogN\og{N + M) ) 



This paper 

o(/W) 


0{aV2\ogN) 


TABLE I 

Support recovery guarantees eor dieeerent sensing matrices. 


note that the random matrix with independent entries is not 
realizable in radar setting but included in the table to provide 
a baseline. 

Uniform recovery guarantees are often formulated in terms 
of satisfying RIP property with high probability, since if 
a measurement matrix satisfies RIP of order 2K such that 
^ 2 K < S ^ \/2 — 1, then all K-sparse vectors are successfully 
recovered with a reconstruction error of an oracle estimator 
that knows the support of the sparse vector or the support of 
K largest elements | [29| . 

Baraniuk et. al. in | [30| have shown that random matrices 
with i.i.d entries from either Gaussian or sub-Gaussian proba¬ 
bility distribution satisfy the RIP condition. For any 5 G [0,1], 
5k ^ ^ K < (^/k) where is dependent on 5 

and the sub-Gaussian norm of the random variables. Although 
these unstructured random matrices have remarkable recovery 
guarantees they do not represent any practical measurement 
scheme, which leads us to consider classical linear time 
invariant (LTI) systems. 

Typically, an active imaging system transmits a signal that 
interacts with a scene of interest and the acquired mea¬ 
surements are used to estimate characteristics of the scene. 
The unknown environment is modeled as an LTI system 
whose transfer function has to be estimated using compressed 
measurements from the data acquisition step. It is assumed 
that there exists a sparse or compressible representation of 
the transfer function in some domain and the goal is to 
solve the sparse estimation problem with the least possible 
measurements. This leads to a structured measurement matrix 
that is either a partial or sub-sampled Toeplitz or circulant 
matrix. The RIP condition of order K for partial Toeplitz 
matrices in the context of channel estimation was established 
by Haupt et. al. in HD- They showed that if the sparsity 
K < then 5k < 5, where depends on 

5. This quadratic scaling of number of measurements with 
respect to sparsity was improved in |[^-p4|. Romberg in 
considered an active imaging system that used waveform 
with a random symmetric frequency spectrum and acquired 
compressed measurements using random sub-sampler or ran¬ 
dom demodulator at the receiver to estimate the sparse scene. 
The resultant system is a randomly sub-sampled circulant 
matrix representing the convolution and compression process. 
It is shown that for a given sparsity level K, the condition 


that 52k < 5 IS satisfied if the number of measurements 
M > aQ5~‘^ min (^K{logN)^,{Klog where > 0 
is a universal constant independent of the size of problem 
and 5. This was extended by Rauhut et. al. in They 
consider a deterministically sampled random waveform in 
time domain with samples following Rademacher distribution, 
which is modeled as a sub-sampled Toeplitz or Circulant 
matrix with entries sampled from Rademacher distribution. 
It was shown that for a given sparsity level 5k ^ 5 
with high probability if the number of measurements M > 
aj mnx {5~^{K log N)^^‘^,5~‘^K{log N log where aj 

is a universal constant. In the subsequent work by Krah- 
mer et. al. in p4| , the relation between sparsity level and 
number of measurements is improved and more general ran¬ 
dom variables are considered such as vectors following sub- 
Gaussian distribution to generate the Toeplitz or Circulant 
matrix. It is shown that, for a given sparsity level K the 
condition 5k ^ 5 is satisfied if the number of measurements 
M > as5~‘^K{logKlogN)‘^, where the constant ag is 
a function of only the sub-Gaussian norm of the random 
variables generating the matrix. We adopt a method similar 
to ED and establish the RIP condition of order K and obtain 
a similar result stating that 5k < 5 if M > a5~‘^K‘^logN, 
where a > 0 is independent of 5. 

The rest of the paper is organized as follows. Section 
gives the mathematical model for the multi-frequency chirp 
model and the statistical model considered for the target. 
Section |nl| states the main result about the measurement 
scheme employed for sparse recovery. Section |IV| contains 
detailed simulation results of the proposed multi-frequency 
chirp waveform. Section |V| contains the detailed proof of the 
main recovery result. We conclude with some future directions 
in Section Ivll 

II. Signal model and problem statement 
A. Multi-frequency chirp model 

We consider a radar sensor with collocated transmitter 
and receiver antennas employing the compressive illumination 
framework proposed in and p4| for estimating the range 
and complex reflectivity of refiectors in the scene. The chirp 
rate of all the transmitted linear frequency modulated (LFM) 
waveform is fixed at where [3 is the bandwidth of each 
transmitted waveform, r is the pulse duration and B = gP is 






























4 


the system bandwidth for > 1. We denote the unambiguous 
time interval as = tmax - imin, where tmax = 
imin = , -Rmax,^min are the maximum and minimum 

range in the area of interest, respectively, while c is the velocity 
of light in vacuum. The whole space of range is discretized 
into grids based on the Radar’s resolution, therefore we get 
N = Btu grids. The interval of frequency from [0, B] is 
divided into N grids such that f{i) = i = 0, • • • ,N — 1, 
which are used as center frequencies for the chirp waveform. 
From the possible N waveform, a subset of size Nc is chosen 
at random for transmission. We simplify this selection model 
by considering independent Bernoulli random variables as 
indicator variables to select LFM waveform such that Nc 
waveform are selected on an average. Let 7 ^ G {0,1} be the 
random variable indicating that f{i) is part of the subset of 
size Nc. It can be seen that 

{ 0 with probability (w.p.) 1 - ^ 

N 

1 w.p. 7 . 

The chosen LFM waveform are then scaled by independent 
random variables such as 

1 ) a sequence of independent and identical complex phase 
with probability density function as ^ 

[0,27r] , 

Ci = 7 *exp(j$j), (4) 


2) a sequence of scaling variables following Rademacher 
distribution given by 


{ —1 w.p. 0.5 
1 w.p. 0.5 

(5) 


We choose the model in which states that the Nc chirp 
waveform are scaled by random signs in our analysis but use 
the model in Q in simulation results. 

The transmitted signal can be written as 


3(i) = 


I iv-l 

7m t. 


Ci exp 




where 0 < f < r. The receiver utilizes stretch processing at 
the same chirp rate as the transmitter and a fixed reference 
frequency fd = fc to demodulate the carrier frequency and 
estimate the round-trip delay. The overall duration of the de¬ 
chirping waveform is + r. The sampling rate employed at 
the receiver is Fg = ^tu. The total number of samples in the 
pulse duration r is M = ptu. The output samples of the stretch 
processor due to the target at different delay bins ^ 

are 


N-l N-1 


y{k) = 


^cjexp -j2» — 




im\ 


N 


X exp ( 27rj ( ^ - ^ ) fc ) x{m) + w{k), 


where k = 0, • • • , M — 1, Wk is measurement noise process 
with 0 mean and variance p = ^ G Z and x (m) is the 


complex scattering coefficient due to a target at the delay bin 
A^. This can be compactly written as 


y = Ax + w, 


( 6 ) 


where y,w G C^, A G and x G C^. The sensing 

matrix A can be represented as a series of deterministic 
matrices corresponding to the response to each of the chirp 
waveform scaled by zero mean random coefficients as shown 


N-l 

A = ^ ^ qH^AD^. 

i=0 


(7) 


The individual components are as follows 


A = 


1 

Vmn; 


[A(0) 


A(7V-1)] 


A(r) = 1 
Y>i = diag 
Hi = diag 


exp(-27rj;^) •• 

• exp 2TTj ^ ^ )) 

1 ^ 

1 exp {-j2TTj^) 

• • • exp j2TT ^ ^ ^ 

)] 

1 exp (i27r]^) • 

•• exp(j27r^7 ^)] 



•) 

( 8 ) 


where i = 0, • • • , A — 1 and r = 0, • • • , A — 1, A G 
are the samples from tones that correspond to each delay bin 
generated as a result of the de-chirping process in case of a 
single chirp system, G is the shift in frequency 

due to the chirp waveform, and G is the phase 

term associated with different delay bins due to the chirp. 
Each column of sensing matrix A can also be represented as 


A(m) = E^FG^c, where 

Em = diag 1 exp (-ji'27r^) • • • exp 

F(r) = [1 exp(27rj^) 


(-j27r^^^) 


exp 5 ) 


Gm = diag 1 exp (-ji'27r^) 


exp 




(9) 


where m = 0, • • • , A — 1, r = 0, • • • , A — 1, G 
represents the tone generated due to target present at delay 
bin, F G are the different chirp center frequencies, 

Gm e is the phase term due to different chirp 

frequencies for a particular delay bin m, and c G is 
the random vector that selects the chirp waveform and scales 
them. A closer inspection of matrix F reveals that each of the 
center frequencies used to shift the chirp waveform is being 
aliased into lower frequency tones as we are sampling at Sub- 
Nyquist rate. We assumed p G Z in order to simplify the 
analysis as we get sub-sampled Discrete Fourier Transform 
(DFT) matrices. We impose an additional condition that p 
should be co-prime with M in order for A frequency tones to 
be uniformly mapped onto M frequency bins, where M < N . 
A simple example of p = M -f 1, which makes p co-prime 
with M, circularly maps the N possible frequencies into M 
bins. 
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B. Target model 

We consider a statistical model similar to Strohmer and 
Friedlander in (H) for the sparse range profile of targets. 
We assume that the targets are located at the N discrete 
locations corresponding to different delay bins. The support 
of the K-sparse range profile is chosen uniformly from all 
possible subsets of size K. The complex amplitude of non¬ 
zero components is assumed to have an arbitrary magnitude 
and uniformly distributed phase in [0, 2tt]. 


C. Problem statement 

Given a sparse scene with targets following the statistical 
model discussed in previous section, and measurement scheme 
in (1^ with M N and sparsity level K N, the goal of 
compressed sensing 0 is to recover the sparse or compress¬ 
ible vector X using minimum number of measurements in y 
constructed using random linear projections. The search for 
sparsest solution can be formulated as an optimization problem 
given below 

min||x||Q, subject to ||Ax — y II 2 < V, (10) 

where rf is the noise variance. This was shown to be NP-hard 
and hence, intractable p5| , and many approximate solutions 
have been found. One particular solution is to use convex 
relaxation technique to modify the objective as norm 
minimization instead of the non-convex ^0 norm given by, 

min||x||^ subject to ||Ax — y Il2<^- (11) 

This approach has been shown to successfully recover sparse 
or compressible vectors given that the sub-matrices 

formed by columns of sensing matrix are well conditioned. 
Our analysis is based on LASSO Q, which is a related method 
that solves the optimization problem in ( pT| ). It has been shown 
in p6| that for an appropriate choice of A and conditions on 
measurement matrix, the support of the solution of the below 
mentioned optimization problem coincides with the support of 
the solution of the intractable problem in 

ininA||x||i + i ||Ax-y|| 2 . (12) 

The goal of our analysis is to show that the measurement 
model given in 0 satisfies conditions on mutual coherence 
given in p6| and to find a bound on the sparsity level of 
range profile, which guarantees successful support recovery of 
almost all sparse signals using LASSO with high probability 
from noisy measurements. In addition, we also provide an 
estimate of the number of measurements required for the 
sensing matrix representing our scheme to satisfy the RIP 
condition. The next section presents our main results of our 
analysis. 


III. Recovery guarantees 

In order to obtain the non-asymptotic recovery guarantee for 
our system employing multiple chirps, we find an estimate of 
the tail bounds of mutual coherence and spectral or operator 
norm of our measurement matrix. Using the estimates, we also 


provide conditions for RIP condition of order K to hold. 

We make use of the Matrix Bernstein inequality given in 
lemma to get a tail bound on the operator norm for the 
measurement matrix given in ([^. 

Lemma 1: Given the measurement matrix model in 0, if 
Nc > ^ log (A -f M), then we can bound the tail probability 
for the operator norm as follows 

P 


HAIL, >2 


Nlog{N + M) 


< 


^N + M 

a2 = — 


a2 — l 


M 


, where 


(13) 


1 + - 


log(Ar+M) 

Ne 


In addition, we also obtain an estimate of the expected value 
of the operator norm of measurement matrix given as 


E 


(llAIL,) < ^J^log{N + M) + 


N 

MNc 

(14) 


The following results about the Euclidean norm of columns 
and mutual coherence are obtained using concentration in¬ 
equalities of quadratic forms of sub-Gaussian random vectors 
given i n p6| , which is extended to the complex domain in 
lemma [12 

Lemma 2: The concentration inequality for the minimum 
of Euclidean norm of any column m of A is given as follows 


^ (imn||A(m )||2 < 1-e) 
™ ( 


< 


exp 


-Md 


\ 




^-1 


(15) 


q* = max |^l,21og 

Lemma 3: If M > (log then there exists constant as > 
0 such that the mutual coherence of our sensing matrix is 
bounded by 


where d > 0 is a universal constant, e G (0,1), and 

N 


pU(a)> 


<T3 

1 — e 


logA\ 

M ) 


< 


+ 4Aexp {-Mde ^), if log A > 

^ + 4 A exp (—Mde^) , otherwise 

(16) 


. N^2 


where d > 0 is a universal constant, e G (0,1) are arbitrary 
constants, and 

2 

. I 

ui = d ' 


U2 = 



d log A, 
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q* = max ( l,21og ( ^ 


e = 


< N ) 


Theorem 1: For a measurement model y = Ax +w, where 
A is defined in 0 such that X is drawn from a K-sparse 
model in complex domain and w ^ CA/'(0, a^J), the following 
conditions guarantee successful support recovery from solving 
with regularizer A = 2cFy/\og N, 

(1 — ei) aiM 


K < Amax = 


log (A) log (A + M) ^ 


M > (logNf ,\ogN > -- 

— V & y 5 to _ y2 /r?*-1) 


Ac > max ( - log(A + M), vN 


0^3 


min \xk\ > 
kei ^/T- 


T^/2^og{N), 


(17) 

(18) 

(19) 

( 20 ) 


with probability P4{1 — pi — P 2 — Ps) for some as > 0, > 

0,ee(0,l), zy^l,d>0isa universal constant independent 
of N,M, where 

exp {-dMf) 

( 1 \“2-l 

= [nTmJ + i-dMe^) 

Ps = 4Aexp (—dMe^) , 

P4 = 1 - 2A-\27rlog (A) + AA-^) - O , 

2 


0(2 — 


ui = d 



g* = max I 1, 2 log 


The proof in section [V| involves direct application of lemma 
and uses the estimates of the spectral norm and the mutual 
coherence of the measurement matrix. 

Theorem 2: For the measurement matrix A given in Q and 
any (5, e G [0,1] such that (5 + e < 1, the RIP condition given 
as < (5 + e is satisfied with probability 1 — ps — Pq if the 
number of measurements M > a6~‘^K‘^ log A, where 
1 


P5 = 


A(^ 3 - 2 ) ’ 

/ 


Pq = 4 A exp 


—d 


N 


M 


Us = a 


N 


a > 0 is a constant independent of A, M. 

We adopt a similar approach as Haupt et. al in (D and 
utilize the estimates of inner-product of columns of sensing 
matrix and norms to obtain a simple bound on the number of 
measurements required to guarantee RIP of order K. 


Discussion 

The support recovery guarantee stated in Theorem is sat¬ 
isfied for almost all K-sparse vectors sampled from the generic 
sparse signal model discussed earlier, i.e. given a measurement 
matrix one could find a K-sparse vector (with arbitrarily small 
probability varying A, M and Ac. ) for which the recovery 
fails. This differs from the worst-case guarantees as well as 
reconstruction error bounds that depend on Restricted Isometry 
Property (RIP) given in Theorem The exponent in probabil¬ 
ity tail bounds for quantities such as mutual coherence /i(A), 
spectral norm || A||^^ of the measurement matrix are controlled 
by number of chirp waveform employed Ac. Specifically, it 
can be seen that the upper bound to the tail probability of 
the above estimates reduce as Ac increases until Ac < ^. 
We also show empirically in section W that the expected 
value of mutual coherence /i(A) reduces as the number of 
chirp waveform increases. Specifically, it converges in mean 
to the mutual coherence of an unstructured random matrix 
G with independent Gaussian entries and thereby converges 
in probability as well. Typically, smaller values of /i(A) are 
desirable for robust recovery as shown by Candes and Plan 
in 1261 as it ensures that the Grammian matrices of the sub¬ 
matrices formed using a subset of columns of sensing matrix 
A are well conditioned as shown by Tropp in p7| . Since, 
the minimum value of the signal has to be above the noise 
floor ( |20| ) for successful recovery, we get a condition on the 
signal to noise ratio SNR for a particular target located at a 
fixed range bin r below which the recovery guarantee does 
not hold, which is given by 

SNRr = > 128 k log AT. 

The authors in GD, 1^ also show that the threshold on 
SNRr scales with Clog A with constant C determining the 
probability of successful recovery.In section |lv| we study the 
effect of SNR on the reconstruction error using simulations. 


IV. Simulation examples 

For our simulation studies, we consider a system with 
a bandwidth B = IGHz from which we choose center 
frequency of each chirp waveform randomly with each chirp 
sweeping a fraction of the system bandwidth ^ We 

note that a wideband multi-tone signal with bandwidth B 
modulating a LFM waveform of bandwidth {3 results in a 
system bandwidth of A + /3 Hz, but we seek to resolve targets 
with range resolution that corresponds to B Hz common to 
all modulated chirps. The fractional bandwidth defined as the 
ratio of ^ represents the under-sampling ratio as the stretch 
processor output is uniformly sub-sampled at that rate. We 
assume that the minumum and minimum range of the area of 
interest are Amin = Om and Amax = 150m , respectively. 
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Fig. 1. Mutual Coherence of structured sensing matrices as a function of 
under sampling ratio. Random Gaussian matrix coherence is provided as a 
baseline. 


The pulse duration r is chosen such that the ratio p = ^ e Z 
is co-prime with the number of samples. We make use of 
the model in to select a subset of Nc chirp waveforms 
and scale with a random phase term to obtain the simulation 
results. The other measurement matrices that we compare the 
performance of our scheme with are 

1) matrix G with i.i.d. complex Gaussian random entries 

sampled from distribution CA/'(0, , 

2) matrix T = which is a partial Toeplitz 

matrix, where Pq is a uniform sub-sampling operator, 


Ti = 


In 

^AT+l tN 


h 

t2 


t2N-l t2N-2 • • • tN. 

ti - CA/'(0,1), i = 1, • • • , 2Ar - 1. and 


We generate 100 realizations of matrices G,A and T, 
in order to observe the effect of under-sampling on mutual 
coherence. It can be seen empirically from figure that the 
mutual coherence of the sensing matrix A representing our 
system converges in mean to the mutual coherence of a sensing 
matrix G as the number of chirps increase. We observe that 
even for small values of z/ = ^ = 0.03, the mean of the 
mutual coherence of our measurement matrix A is quite close 
to the mean of the mutual coherence of matrix G. In addition, 
we also evaluate the coherence of Toeplitz matrix T, which 
is representative of active-imaging schemes using stochastic 
waveform that are modeled as linear time invariant systems. 

Next, we consider the recovery performance of our mea¬ 
surement system used in conjunction with Basis pursuit de¬ 
noise using SPGLl solver developed by Van den Berg and 
Friedlander in p^ , to estimate the unknown target 

locations and their amplitudes in the area of interest. For each 


realization of the measurement matrix we generate multiple 
samples of target range profile with specified sparsity level 
and scattering coefficient is sampled at specified locations are 
sampled from a complex Gaussian distribution. The overall 
target range profile is normalized to get a fixed SNR. We 
consider a function of mean square error as a performance 
measure, specifically a thresholding function for a fixed SNR 
of 25dB. We consider a threshold of 1% on the mean squared 
error and vary both the number of targets in the scene, and 
the bandwidth of the chirp waveform, which in turn infiuence 
the sampling rate at receiver. Again it is clear from figure 
that the performance becomes similar to that of the random 
Gaussian sensing matrix G as the number of chirps increase. 
We also observe that the recovery performance of Toeplitz 
matrix T is similar to our system at lower values of fractional 
bandwidth I <0.7 but performs better when I > 0.7. 

Next, we see the infiuence of noise on sparse target recovery 
using Basis pursuit de-noise employing our measurement 
scheme. We fix fractional bandwidth ^ = 0.4 and vary 
the noise variance as well as the number of targets in the 
scene. In Figure the intensity of the image represents the 
mean square error in dB scale. Figure shows that as the 
number of chirps increase, the performance achieved by our 
scheme in terms of the mean square error approaches the 
mean square error achieved by the random Gaussian matrix G. 
The reconstruction error for Toeplitz matrix T is marginally 
better at lower values of SNR « l^dB compared to our 
measurement scheme. 


V. Proofs 


Proof of lemma Uk The result can be obtaine d by d irect 


application of lemma 


using the value of L = 


N 

MN^ 


and 

z^(A) = ^ obtained from lemmaj^and lemma|^ respectively. 
The upper bound on the tail probability is given as 


(||Atp>t) <(7V + M)exp 


-t‘^/2 


N t I TV 
^ M , 


MAT, 


By plugging in t = we get the result in ( p3] ). 

For the tail probability to decay, w e require t hat 0^2 > 1. 

This gives us the condition that <0.5, which 

implies Nc > ^ log {N + M). Similarly, using the estimates 
L, and o{A) we can bound the expected value of the operator 
norm of A as given in lemma ■ 

Proof of lemma The norm of a column m of the 
sensing matrix A can be written as follows 

||A(m )||2 = c*Bc, 

where B = G^F*E^E^FG^, and c G is a sequence of 
random variables that selects and scales a subset of the chirp 
waveforms. It can be verified that the diagonal elements of 
matrix B are as follows 

Bi^i = = !,■■■ ,N. 

































(a) A with Nc = 5 Chirps (b) A with Nc = 25 Chirps (c) Gaussian sensing matrix G. (d) Toeplitz sensing matrix T. 


Fig. 2. Phase transition curve for the multi-frequency measurement scheme as number of chirps Nc increases and other related sensing schemes. Image 
intensity is either 1 if reconstruction error is below 1% and 0 otherwise. We compare the performance of our multi-frequency scheme with other random 
structured and unstructured matrices as the sparsity level and fractional bandwidth is varied at Signal to Noise ratio of 25 dB. 



Signal to Noise Ratio in dB 



Signal to Noise Ratio in dB 



Signal to Noise Ratio in dB 



(a) A with Nc = 5 Chirps 


(b) A with Nc = 25 Chirps 


(c) Gaussian sensing matrix G (d) Teoplitz sensing matrix T 


Fig. 3. The intensity values of the images represent the reconstruction error in the log scale. Each image from left to right, represents the recovery error for 
our measurement scheme as number of chirps Nc is increased and other measurement schemes as a function of SNR and sparsity level at a fixed fractional 
bandwidth ^ = 0.4. 


Since the random variables Ci are independent and E (c^) = 0, 
the off-diagonal terms vanish and we get 


E 


i=l 

1 ^ 

= ivg'’ 

E(||A(m)||2) =1. 



- 1 


) 


> f < 


nA^) 1 r^i 

N Nr-a* I M \ 


K 




M rN- 


_r2v 1 

N NcQ* I M I 


< 1 




< 


4 exp — 


Nciq*) 1 rN_i 
N Nr-a* I M I 


mm 


NcQ* \ M 


iVe(^). M rN- 


n 


* \ M \ 


NcQ* I M 




e (0,1]. Using the 


Let t = e for any e < 
approximation ~ we get the required result. ■ 
Proof of lemma We can express the inner-product 
between any two columns mi and m 2 of sensing matrix as 

(A(mi), A(m 2 )) = c^Bc, 

where B = G^^F*E^^E^ 2 FG^ 2 - The diagonal terms of 
the matrix B is given as 

'mi — m 2 \ 1 


= D 


M 


Using results from lemma and lemmaalong with the 
result on sub-Gaussian norm from lemma Mil we have 


N 


Dm 


mi — m 2 

N 


Nc 

M-l 


exp 


27rjf(mi — m2){i — 1) 

N 




t=0 


N 


where i = !,••• ,A. By using the fact that q are zero 
mean independent random variables, we obtain the following 
expression for E (c^Bc) 


( 21 ) 


E(c^Bc) =^E(cf)Bi^ 


= Dm 

= 0 . 


mi — m 2 

N 


N 


^exp 

i=l 


2TTj{mi - m2){i - 1 ) 

N 


The concentration inequality for the minimum value of norm 
of any column can be written as for any t G [ 0 , 1 ], 


P ^min ||A(m)||2 > I — > 1 — NP ^||A(m)||2 < 1 — , 

P (mn||A(m )||2 < 1 - i) < NP(\\A{m)\\l < 1 -. 


Using the results from lemma and lemma along with 
the sub-Gaussian norm result from lemma and making the 
approximation ~ we see that there 3 0^3 > 0 such 
that 


P |(A(mi), A(m 2 ))| > as 


log A 
M 


< 
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4exp I I , 


Nc\^ 
N 


(f-1) 


M 


Where h{N) = min ^ y lii^Ar 

M > (logA^)3 and \§]- we get 


. Using the fact that 


P |(A(mi), A(m 2 ))| > as 


log 

M 


< 


iv^r if logfV > 
Otherwise, 


4 

. N'^2 


( 22 ) 


where 


4^1 = (i 


q c ^3 


Ac VC, 

AT 


(^-1) 


^2 = — . , dlogA^. 

A^c \ q* / 

N 


/i(A) = max 


mi,m2 ||A(mi)|| ||A(m2)|| 

< max |(A(mi), A(m 2 ))| max- 


1 


mi ,m2 


A(r 


^U(A)>- 


0^3 /log(A)\ 


M I 


<p( max |(A(toi), A(to 2 ))| > I 

\ mi,m 2 \ M I 

+ P ^min ||A(m )||2 < 1 — 

|{A(mi), A(to 2 ))| > 

+ P ^min ||A(m )||2 < 1 — 

Using and ( | 22 | ) in the above expression, we get the result 
in ( p^ . ■ 

Proof of Theorem^ Using M > log(A)^ and log A > 

” I 26 



iVcU 2 / *- 1 ) 0^3 coherence condition given in 

(iv" 

13 satisfied with high probability as shown below 
1 


liiA) = 0 


W.p. Pi > 1 


log N 




+ 4Aexp (—dMe^) , 


where 0^3 > 0 is a constant independent of N and M, e > 0, > 


ui = d 


r^3 


Ac Vq 
N 




_ 

tail bound in ( p^ depends on ^ as the function 
is an increasing in Nc G [ 0 , l/exp(l)] but decreases in 
Nc G [l/exp(l), 1]. In addition to this, we also believe that 
0^3 decreases as Nc increases linearly with N, which is also 
verified numerically in section This leads to part of the 
condition in ( p^ given as ^ ^ 1 . 

The measurement matrix in our analysis does not have unit 
norm columns and in order to apply lemma we normalize 
the columns. We follow the approach similar to GD Let D G 
j^ATxAr (diagonal matrix with diagonal entries corresponding to 
the norm of the column of A given by 

A,4 = ||A(z)|| 2 . 

The measurement model can be modified as 
y = Az + Av, 

where A = AD“^ and z = Dx. Next, we obtain the prob¬ 
ability tail bound for the operator norm of the measurement 
matrix A with £2 normalized columns. Using lemma we 
have Ve > 0, e G (0,1), independent of N and M, 


The concentration inequality for the coherence of matrix A 
can be obtained by using the following inequality 

|(A(mi), A(m2))| 


> 


\/l — I 


A log {N PM) 
M 


< P 


l|A|L^ > 2 


A log {N PM) 
M 


P P (imiii. Dcci A c j 

\ m J 


< 


N PM 


ai — 1 


+ 4Aexp (—dMe^) , 


where 


2 / log(Ar+M) 


Qfl = 

f + iv n;_ 

A, >Log(A + M), 


(t) 


A-1 


Therefore, 


< 


op 


a/T^ 


A log [N PM) 


w.p. P 2 > 1 - 


1 


M 

ai — 1 


(24) 


m exp {-dMe^) . 


,N + M^ 

This gives us the condition in ( fT^ . Using in iSZi we 


(23) 


and g* = max ^1,2log This establishes the condition 

in ( p^ . We also note that the exponent in the probability 


obtain the condition in ( p^ . 

Next, we establish that the measurement matrix does not 
reduce the absolute value of non-zero entries of the sparse 
vector X below the noise level. 

P (jmnDi^i \xi\ < 8 ( 7 a/ 2 log A^ < NP [Di^i < y/l — e) 

< 4Aexp(-dMe^) . 

Therefore, we have 

min l^il > SapApN (25) 
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w.p. P 3 > I — 4A^exp (—. 

We define the following events associated with a realization 
of measurement matrix A 

E,: 

op Amax loS N 

53 : min \zi\ > ^(7^/2\ogN. 

i 

54 : successful support recovery for a fixed sensing matrix. 


Let S be the event that the sampled measurement matrix 
satisfies the conditions required for successful recovery and 
recovers a K-sparse vector x selected from the target model. 
This implies 

p(s)>P(S4|SinS2nS3)x 

(l-P(SJ)-P(S-)-P(S-)). (26) 

Using result from Lemma for P (S4 I Si n S2 n S3),([24|), 
( [^ and ( |^ in ( [^ , we get the desired recovery guarantee. 


Proof of Theorem ^ Using 5k < (5 + e in (|^, it 
can be deduced that the eigen values of ApAr,Vr are 
G [1 — (5 — e, 1 + (5 + e] such that card (L) < K. This 

can be translated into a condition on the elements of the 
Gramian matrix A* A using lemma given as the follow- 
ing events Sg : {|(A(mi), A(m2))| < 4} Sg : 

mi/m2 


n™=i{|l|A(m)||^-l|<e}. 

Let S be the event denoting that the RIP condition of order 
K is satisfied. Then we have 

P(S4 <1-^P (KA(mi), A(m2))| > ^ 

-AP(|||A(m)||2-l| >e) (27) 

Using the result from ( [21] ) and using the fact that we can 


choose Nr such that e < 




, we get 


( I|A(to )||2 - 1 > e) < 4exp 


—d 


(^-1) 


N 



Similarly, using results from lemma [^ we can obtain 


P ( |(A(mi), A(m2))| > ] < 


4 exp 


V 


-"4 


N 



(29) 


Using ( [^ , ( |^ and the condition that M > aS ‘^K‘^\ogN 
in we get the desired results. ■ 


VI. Conclusion 

In this work we have shown that structured compressed 
sensing matrices with random components can be realized 
in radar (delay estimation) setting using an LFM waveform 
modulated by a random sparse multi-tone signal in transmit 
and a simple traditional analog receiver structure. We provide 
recovery guarantees for the proposed compressive sensor com¬ 
parable to random Toeplitz/Circular matrices with much larger 
number of random elements. The proposed scheme is well 
matched to practical implementation utilizing small number of 
random parameters and uniform sampling ADCs on receive. 

A potential direction for future research is to investigate the 
effectiveness of multi-chirp waveforms in the multiple output 
multiple input (MIMO) setting with multiple transmit and 
receive antenna elements for estimating support of targets in 
angle and range domain (m We note that in this setting each 
transmitter can use a single chirp with a random frequency 
offset and superposition is achieved at each receiver as the 
waveforms reflected from the scene is naturally summed at 
each receiver. 


Appendix 

Some Useful Lemmas 

Lemma 4: Given a complex matrix M G with Ci = 

Mi^i, and Ri = \Mij\ then we have 

A, G U^=iD(Q, Vi = 1 ,... , n, ( 30 ) 


where D(c, r) is a disc with center c and radius r, and are 
the eigen-values of M. 

We restate the theorem given in | [T3| , which was extended to 
the complex setting in 1^ , that give the conditions on the 
measurement matrix for successful recovery when using £i 
penalized optimization methods. 

Lemma 5: For a measurement model y = Ax + w, where 
A G has unit £2 norm columns, x is drawn from a 

K-sparse model in complex domain and Wi ^ Cff{0,cF‘^), if 
the following conditions are satisfied 


Ai(A) < 


ao 

log A’ 


(31) 


where ckq > 0 is a constant independent of the dimensions of 
the problem; also, if 


K < (A) 


aiN 

l|A||^,log(A)’ 


(32) 


for some o^i > 0, and 


min|a:;fe| > 8a\/2 log (A), (33) 

kei 

then the solution x of Gl has the same support as the 
unknown sparse vector x and relative error is bounded as 
shown below 


supp{^) = supp{x.) (34) 

l|x-x||^ aVm 

I|x|l2 ” Ily|l2 

> l-2A-427rlog(A) +AA-i)-0(A“^'°s^). (35) 


























11 


Lemma 6 (Matrix Bernstein inequality Let be a 

sequence of i.i.d. random matrices. For a random matrix 
expressed as A = A^ we have 

P (llA^p > < {di + di) exp ^ +/(A) ) ’ 

l|Aitp<i:,Vz = i,--- ,z) 

z/ (A) = max (E (AA*), E (A*A)), 

where A^ G The expected value of the operator norm 

of A is bounded by 

E (||A||^p) < ^2v (A) log(di +^2) + +^2) ^ 

(37) 


Lemma 7: Given a matrix A, the operator norm of the 
matrix A = AD-i is bounded by the following inequality 


mim Dj 


> 


> 


op 


max. Dj 


(38) 


where D is a diagonal matrix with positive diagonal elements. 
Proof: For any vector v such that ||v||2 = 1, we have 


Av = 

^ l|Av||. 

where u(v) = 


D ^vll = Au(v) IID ^vI 


= l|Au(v)||, 


l2 
|D-^ 


D-"v 

l|D-iv|L 


-. Since ||v|| 
Euclidean norm of as follows 

1 


1, we bound the 


min^ Dj 

Using this we get, 


> D- 


I2 - max, A 


l|Au(v)||2 > ||Av|| > 

min, A,i ~ II2 — max, 

By taking supremum over v in the space of unit norm vectors 
we obtain the desired result. ■ 

Lemma 8: For the matrix c,H,AD, G in ([^, we 

have 


gH, AD, 


< 
op — 



(39) 


Vi = I,-- - 

Proof: By sub-multiplicativity property of the operator 
norm we have, 


||c,H,AD,||^^ < Id ||H,tJ|A||^J|D,t^. 

Since H, and D, are diagonal matrices with complex expo¬ 
nential entries, it can be shown that ||H,||^^ = ||D,||^^ = 1. 
Also, we assume that c, are sub-Gaussian random variables 
with |c,| < 1. In order to find the operator norm of A, 
we define G = A A* G since it is full rank and 

II “ 112 

A = IIGII The entries of matrix G are as follows 

II 11 op " ''Op 


G{k,k) 

G{k,l) 


N 

mW, 

N 1 
MNcN 


N-l 

^exp 

m=0 




Dn ( 1 = O5 


MNr 


N 


\/k, I = 0, • • • , M — 1, such that k ^ I, and (^) = 
^ exp (27ri ^m) is the discrete Dirichlet Kernel. The 

second term is zero because the discrete Dirichlet kernel is 
being evaluated at it’s zeros, which are the Fourier frequency 
bins {^),n eZ. Therefore, 

N 


G = 


MK 


I. 


This implies that the ||G|| = and leads to the result in 

([33. ' ■ 

Lemma 9: For the matrices P, = c,H,AD, G 
given by ([^, we have 


N 


i=l 

N 

^E(PiP*) 


N 

< —, 
- M’ 




N 

M' 


(40) 


(41) 


Proof: First, we compute the norm of P,P* as it is a full- 
rank matrix using D,D* = I, H,H* = I and A A* = 
we have 


P,P* = c,c*H,AD,D*A*H* 
. N ^ 


' ^ MNc ' 

Using the probabilistic model for c, given in ([^, we get 


E(c.c*) = f. 


We have 


N 




5:e(p.p-)=j^l 




Applying the operator norm yields the result in ( |4T] ). Sim¬ 
ilarly, using the sub-additivity of the operator norm and 
||H,AD,D*A*H*|| = ||D*A*H*H,AD,|| we get 

11 t I I ^ Wop 11 7 2 * MI op ® 


N 




2 = 1 


op 


N 

^E(c*Ci)D*A*H*HiAA 
2 = 1 
N 

<^|Edc*)|||H,ADiD*A*H*| 
2 = 1 
_ N 

~ M' 


op 


then we have 


Lemma 10: Let B = G^F*E;^E^FG^, and B = 


I|BL, < 


B < 

I 11 op — 


N 

M 

'N 

M 


A’ 


N 

M 

N' 

M 


Nc 

\/M 


(42) 

(43) 


where [re] = 2; G Z such that z > x, Vre G M. 
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Proof: We note that RankCB) = RankCB) = M, and 
we can obtain a bound on the Frobenius norm of the matrices 
as shown below 


In order to find the bound on the operator norm we see that 



Since and are diagonal matrices with complex 

exponentials along the principal diagonal, it can be shown that 
lIGmllop = IIEmIlop = 1 - In order to estimate the ||F||^p, we 
see that 



■ Vi,i ■ 


F*F = — 

V2,l • 


Nc 


-1 

cT 

> 


Since p is co-prime with M, we observe that the N possible 
frequency tones circularly get mapped onto M possible aliased 
sinusoids. Therefore, if neither i,j 7 ^ 




If either i = \%] or j = then Vij e is a partial 

identity matrix. It can easily be verified that ||F||^^ = 

Using this result along-with the estimate on the bound for the 
Frobenius norm, we get the desired results. ■ 

Lemma 11: For the sub-Gaussian random variables de¬ 
scribed in 0 and (|^, the Sub-Gaussian norm ED is as 
follows 


IlCilU, = supEdCiH^ -L 


g>i 


Nr. 


q* = max U, 21 og 


(44) 


Proof: For the probability models given in ^ and 
we have 


0 with probability 1 — ^ 

1 w.p. 




1 


The solution to this optimization problem can be found by tak¬ 
ing the logarithm and solving the unconstrained optimization 
problem which is given as 



In order to satisfy the constraint, the solution is lower bounded 
by 1 . ■ 

Lemma 12: Given a zero mean real random vector c com¬ 
posed of independent and sub-Gaussian random variables 
Ci,i = I,-- - ,N such that ||ci||^^ < K, i = 1, - ■■ ,N, wc 
have 

Pr (|c^Bc-E(c^Bc)| > t) 

where B e , c € for some absolute constant d> 0 

and Vf > 0. 

Proof: Let B = B/j + iB/^, where B/j, B/^ G . 

Therefore, using the Hanson-Wright inequality for real ma¬ 
trices given in and the fact that ||Bij||^p < ||B||^p, 
WBi^Wop < IIBII^and WB^Wp < ||B||^, ||B,„||^ < ||B||^ 
we have 

Pr (|c^Bc - E (c^Bc) | > f) < 

Pr ^|c^Bi{C - E (c^Bijc) | > + 

Pr ^|c^B/„c - E (c^Bimc) \ > ■ 

This gives us the inequality in ( [45] ). ■ 

Acknowledgment 

This research was partially supported by Army Research 
Office grant W911NF-11-1-0391 and NSF Grant IIS-1231577. 

Reeerences 

[1] M. A. Richards, Fundamentals of Radar Signal Processing. New York, 
NY, USA: McGraw-Hill, 2014. 

[2] R. Middleton, “Dechirp-on-receive linearly frequency modulated radar 
as a matched-filter detector,” IEEE Transactions on Aerospace and 
Electronic Systems, voL 48, no. 3, pp. 2716-2718, July 2012. 

[3] D. L. Donoho, “Compressed sensing,” IEEE Transactions on Informa¬ 
tion Theory, vol. 52, no. 4, pp. 1289-1306, Apr. 2006. 

[4] E. Candes, J. Romberg, and T. Tao, “Robust uncertainty principles: exact 
signal reconstruction from highly incomplete frequency information,” 
IEEE Transactions on Information Theory, vol. 52, no. 2, pp. 489-509, 
Eeb. 2006. 

[5] E. Candes and T. Tao, “The dantzig selector: Statistical estimation when 
p is much larger than n,” The Annals of Statistics, vol. 35, no. 6, pp. 
2313-2351, Dec. 2007. 

[6] S. S. Chen, D. L. Donoho, and M. A. Saunders, “Atomic decomposition 
by basis pursuit,” SIAM Rev., vol. 43, no. 1, pp. 129-159, Jan. 2001. 

[7] R. Tibshirani, “Regression shrinkage and selection via the lasso,” Jour¬ 
nal of the Royal Statistical Society: Series B, vol. 73, no. 3, pp. 273-282, 
1996. 

[8] J. Tropp and A. Gilbert, “Signal recovery from random measurements 
via orthogonal matching pursuit,” IEEE Transactions on Information 
Theory, vol. 53, no. 12, pp. 4655-4666, Dec. 2007. 

[9] D. Needell and J. Tropp, “Cosamp: Iterative signal recovery from incom¬ 
plete and inaccurate samples,” Applied and Computational Harmonic 
Analysis, vol. 26, no. 3, pp. 301 - 321, 2009. 

[10] L. Potter, E. Ertin, J. Parker, and M. Cetin, “Sparsity and compressed 
sensing in radar imaging,” Proceedings of the IEEE, vol. 98, no. 6, pp. 
1006-1020, June 2010. 

[11] K. Gedalyahu and Y. Eldar, “Time-delay estimation from low-rate 
samples: A union of subspaces approach,” IEEE Transactions on Signal 
Processing,, vol. 58, no. 6, pp. 3017-3031, June 2010. 

[12] M. Herman and T. Strohmer, “High-resolution radar via compressed 
sensing,” IEEE Transactions on Signal Processing, vol. 57, no. 6, pp. 
2275-2284, June 2009. 






13 


[13] T. Strohmer and B. Friedlander, “Analysis of sparse mimo radar,” 
Applied and Computational Harmonic Analysis, vol. 37, no. 3, pp. 361 
- 388, 2014. 

[14] T. Strohmer and H. Wang, “Accurate imaging of moving targets via 
random sensor arrays and kerdock codes,” Inverse Problems, vol. 29, 
no. 8, p. 085001, 2013. 

[15] M. Hgel, H. Rauhut, and T. Strohmer, “Remote sensing via ii- 
minimization,” Foundations of Computational Mathematics, vol. 14, 
no. 1, pp. 115-150, 2014. 

[16] M. Rossi, A. Haimovich, and Y. Eldar, “Spatial compressive sensing in 
mimo radar with random arrays,” in Information Sciences and Systems 
(CISS), 2012 46th Annual Conference on. Mar. 2012, pp. 1-6. 

[17] M. Shastry, R. Narayanan, and M. Rangaswamy, “Sparsity-based signal 
processing for noise radar imaging,” IEEE Transactions on Aerospace 
and Electronic Systems, vol. 51, no. 1, pp. 314-325, Jan. 2015. 

[18] E. Baransky, G. Itzhak, N. Wagner, I. Shmuel, E. Shoshan, and Y. Eldar, 
“Sub-nyquist radar prototype: Hardware and algorithm,” IEEE Transac¬ 
tions on Aerospace and Electronic Systems,, vol. 50, no. 2, pp. 809-822, 
Apr. 2014. 

[19] J. Yoo, C. Turnes, E. Nakamura, C. Le, S. Becker, E. Sovero, M. Wakin, 
M. Grant, J. Romberg, A. Emami-Neyestanak, and E. Candes, “A 
compressed sensing parameter extraction platform for radar pulse signal 
acquisition,” IEEE Journal on Emerging and Selected Topics in Circuits 
and Systems, vol. 2, no. 3, pp. 626-638, Sep. 2012. 

[20] J. Tropp, J. Laska, M. Duarte, J. Romberg, and R. Baraniuk, “Beyond 
nyquist: Efficient sampling of sparse bandlimited signals,” IEEE Trans¬ 
actions on Information Theory, vol. 56, no. 1, pp. 520-544, Jan. 2010. 

[21] M. Duarte and R. Baraniuk, “Spectral compressive sensing,” Applied 
and Computational Harmonic Analysis, vol. 35, no. 1, pp. Ill - 129, 
2013. 

[22] A. Gilbert, M. Strauss, and J. Tropp, “A tutorial on fast fourier 
sampling,” IEEE Signal Processing Magazine, vol. 25, no. 2, pp. 57-66, 
Mar. 2008. 

[23] E. Ertin, “Frequency diverse waveforms for compressive radar sensing,” 
in Waveform Diversity and Design Conference (WDD), 2010 Interna¬ 
tional, Aug. 2010, pp. 000216-000219. 

[24] E. Ertin, L. Potter, and R. Moses, “Sparse target recovery performance of 
multi-frequency chirp waveforms,” in 19th European Signal Processing 
Conference, 2011, Aug 2011, pp. 446-450. 

[25] T. T. Cai and T. Jiang, “Limiting laws of coherence of random matrices 
with applications to testing covariance structure and construction of 
compressed sensing matrices,” The Annals of Statistics, vol. 39, no. 3, 
pp. 1496-1525, June 2011. 


[26] E. Candes and Y. Plan, “Near-ideal model selection by 1 minimization,” 
The Annals of Statistics, vol. 37, no. 5A, pp. 2145-2177, Oct. 2009. 

[27] K. Davidson and S. Szarek, “Local operator theory, random matrices and 
banach spaces,” Handbook on the Geometry of Banach spaces, vol. 1, 
pp. 317-366, 2001. 

[28] W. Bajwa, “Geometry of random toeplitz-block sensing matrices: bounds 
and implications for sparse signal processing,” Proc. SPIE, vol. 8365, 
pp. 836505-836505-7, 2012. 

[29] E. Candes, “The restricted isometry property and its implications for 
compressed sensing,” Comptes Rendus Mathematique, vol. 346, no. 9- 
10, pp. 589-592, May 2008. 

[30] R. Baraniuk, M. Davenport, R. A. DeVore, and M. B. Wakin, “A 
simple proof of the restricted isometry property for random matrices,” 
Constructive Approximation, vol. 28, no. 3, pp. 253-263, Dec. 2008. 

[31] J. Haupt, W. Bajwa, G. Raz, and R. Nowak, “Toeplitz compressed 
sensing matrices with applications to sparse channel estimation,” IEEE 
Transactions on Information Theory, vol. 56, no. 11, pp. 5862-5875, 
Nov. 2010. 

[32] J. Romberg, “Compressive sensing by random convolution,” SIAM J. 
Img. Set, vol. 2, no. 4, pp. 1098-1128, Nov. 2009. 

[33] H. Rauhut, J. Romberg, and J. A. Tropp, “Restricted isometries for par¬ 
tial random circulant matrices,” Applied and Computational Harmonic 
Analysis, vol. 32, no. 2, pp. 242 - 254, 2012. 

[34] F. Krahmer, S. Mendelson, and H. Rauhut, “Suprema of chaos processes 
and the restricted isometry property,” Communications on Pure and 
Applied Mathematics, vol. 61, no. 11, pp. 1877-1904, 2014. 

[35] S. Foucart and H. Rauhut, A Mathematical Introduction to Compressive 
Sensing. Birkhauser Basel, 2013. 

[36] M. Rudelson and R. Vershynin, “Hanson-wright inequality and sub- 
gaussian concentration,” Electron. Commun. Probab., vol. 18, pp. no. 
82, 1-9, 2013. 

[37] J. Tropp, “Norms of random submatrices and sparse approximation,” 
Comptes Rendus Mathematique, vol. 346, no. 2324, pp. 1271 - 1274, 
2008. 

[38] E. van den Berg and M. R Eriedlander, “SPGLl: A solver for large-scale 
sparse reconstruction,” June 2007, http://www.cs.ubc.ca/labs/scl/spgll. 

[39] -, “Probing the pareto frontier for basis pursuit solutions,” SIAM 

Journal on Scientific Computing, vol. 31, no. 2, pp. 890-912, 2008. 

[40] J. Tropp, “An Introduction to Matrix Concentration Inequalities,” ArXiv 
e-prints, Jan. 2015. 

[41] R. Vershynin, Introduction to the non-asymptotic analysis of random 
matrices -Compressed sensing, 210-268. Cambridge, USA: Cambridge 
Univ. Press, 2012. 



